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Abstract 



■ The parity-non-conserving (pnc) asymmetry in n+p — ► cZ+7 at thermal energies 

VjD . has recently been calculated using effective field-theory methods. A comparison of 

this calculation with much more elaborate calculations performed in the 70's is 
made. This allows one to assess the validity of this new approach as presently used. 
It is found to overshoot the almost exact calculations by a factor close to 2 for the 
contribution involving the 3 Si component of both the initial and final states. This 
is much larger than anticipated by the authors. This discrepancy is analyzed and 
found to originate from the over-simplified description of the deuteron and capture 
states which underlies the new approach. The claim that earlier determinations of 
the sign would be in error is also examined. It is found that the sign discrepancy is 
most probably due, instead, to the fact that the pion-nucleon interaction referred to 
by the authors corresponds to a parity-non-conserving potential with a sign opposite 
to what is currently used. Some estimates and constraints relative to the pnc 7rNN 
coupling, h^, which the above asymmetry is dependent on, are reviewed. Further 
details are given in an Appendix. 



1 Introduction 

At the time where a new measurement of the parity-non-conserving (pnc) asymmetry 
in radiative neutron-proton capture at low energy is in preparation at LANSCE Q, in- 
tending to improve upon a previous one at ILL 0, it is appropriate to wonder about 
the validity of theoretical estimates of the effect. This one is expected to be dominated 
by the contribution of the pnc pion-exchange force and its measurement should therefore 
provide unique information on the pnc 7rNN coupling, h\, which it strongly depends on. 
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This coupling constant is an important ingredient of pnc NN forces and its determination 
represents a major goal in the field (see |3j for a review). This program supposes that 
theoretical uncertainties related to the description of the deuteron and np capture states 
are small enough so that the coupling can be reliably extracted from the measurement. 
By examining studies performed in the 70's [§J |, || [7j and later ||, some values with 
different signs and magnitudes can be found. Yet, it was concluded that the estimates of 
the asymmetry were rather independent of the strong interaction model |§ but not much 
detail was published on the relevance of some ingredients (see below). When correct ones 
are used, the relation between the pnc asymmetry to be measured, Ay, and the above 
coupling, hi, is essentially given by: 

Ay = -0.11/4. 

From the comparison of results obtained with reasonable strong interaction models, an 
uncertainty of about 5% could be ascribed to the above value. This one is slightly larger 
than the first estimate by Danilov [[TIJ, Ay = — 0.08 /i*, which was obtained on the basis 
of an approach directly relying on nucleon-nucleon phase shifts, but neglecting the effect 
of the tensor force. 

Recently, a new calculation based on a different approach, namely an effective field 
theory [11|, came with a larger value, Ay = +0.17 h\, and furthermore an opposite sign. 



An uncertainty of about 30% was ascribed to this value by the authors. This one could 
thus be compatible with the previous one, but it could also be twice as large. Apart from 
the fact that the discrepancy invalidates at first sight the argument that the coupling, 
hi, could be determined unambiguously from the measurement of Ay, it is important 
to determine how reliable is the new calculation and how it should be corrected when 
improvements are made. Notice that values of the asymmetry with a size comparable to 
the above one may be found in the literature but they correspond to a singlet scattering 
length inappropriate for the neutron-proton system. 

The sign of the effect is another issue which deserves attention. It is the first information 
that one can hope to extract from an experiment, prior to the size itself. When calculations 
of the asymmetry Ay were performed in the early 70's, it was believed that the sign of the 
pnc coupling, hi, could not be determined theoretically, partly explaining the existence 
in the literature of relations of Ay to hi with different signs ||. Since a theoretical 
determination of this sign has been made relatively to the parity conserving one, g^NN, 
(this is the information that matters) |T2]] , it is essential to establish what is the correct 
sign for the relation of Ay to hi. 

It this paper, we intend to compare the recent calculation to earlier ones, which, because 
they were quite sophisticated, can be used as a benchmark for testing the applicability 
of the effective field-theory methods to the process under consideration. At the same 
time, we discuss the points that are sources of discrepancy between calculations. After 
reminding a few ingredients required for the calculation of the pnc asymmetry of interest 
here, we show that the expression obtained in the new approach at the leading order 
identifies to a standard potential based calculation employing zero-range nucleon-nucleon 
forces (Sect. 2). In the third section, we discuss the size of the effect in relation with a 
better description of the nucleon-nucleon states. This is partly done by using analytical 
calculations, which has the advantage to show how the result changes while improvements 
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are made. Section 4 is devoted to the sign of the effect. Details are given, allowing 
one to check its derivation. We, in particular, discuss the independence of this sign 
with respect to conventions that appear at intermediate steps of the calculation. In the 
fifth section, we present a short critical review of what is known, both theoretically and 
phenomenologically, on the elementary pnc 7rNN coupling constant, h\. More details 
about this coupling are given in an Appendix. 



2 Main ingredients entering the calculation of h\ 



Due to parity-non-conservation in nucleon-nucleon forces, the emission of photons in the 
thermal neutron-proton capture can evidence an asymmetry with respect to the neutron 
polarization. This one, Ay, is defined by the angular dependence of the capture cross 
section on the angle, 9, between the direction of the photon emission and the neutron 
polarization, W{9) oc (1 + Ay cos 0). A simplified expression, which accounts for the 
dominant contributions of interest here, is given by: 



In this formula, the factor 2 is a current one appearing in expressions for an asymmetry. 
The sign is well defined but supposes that appropriate conventions are used to calculate 
the matrix elements, El and Ml. This will be detailed in the fourth section together 
with the spin-isospin algebra relative to the corresponding operators as well as to the pnc 
pion-exchange potential. The calculation of the two matrix elements are performed from 
the standard expressions for the electric and magnetic, A J = 1, transition operators. 
For the first one, we rely on the Siegert theorem, which allows one to sum up various 
contributions implied by gauge invariance, including those pertaining to the description 
of the pnc interaction. The sign difference which depends on whether the photon is 
absorbed or emitted is taken care of. For simplicity, factors that cancel out in the ratio 
given by Eq. (1) are ignored. 

The matrix element Ml of interest here involves the overlap of the wave functions of 
the deuteron state and the 1 5'o scattering state, which is the dominant one and provides 
most of the contribution to the capture cross section. It is given by: 



The matrix element, El, involves a transition between the deuteron state and the 3 Si 
scattering state. Its value, which is zero if parity is conserved, depends on the odd-parity 
component admixed to either state by the pnc nucleon-nucleon force. Capture from 
the 3 Si state is not large, but it appears here because a state with a polarized neutron 
necessarily results from its coherent superposition with a 1 Sq state. The expression of the 
matrix element is given by: 




(2) 




(3) 
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where the factor | comes from integrating over angles the angular dependence due to 
both the dipole operator and the parity admixed component with the result 

dil , , 1 



/ — r V = - r z 5 3 . 
J An 3 



The minus sign in Eq. (3) comes from taking the complex conjugate of the complete 
deuteron wave function whose parity admixture contains the extra imaginary number, 
i (which combines with a i in the electric dipole operator to produce a real number). 
For simplicity, the admixture of a 3 Di component to the 3 Si one has been omitted in the 
above equation (as it was in but it is accounted for in the most elaborate calculations 
whose results are reminded below. For the calculation of interest here, its effect can be 
accounted for by replacing the current 3 -wave component, u(r), by the combination 
involving the D-wave component, u(r) + -^w(r). As to the quantities, Vrf( r ) an d 
they respectively describe the parity admixed components to the deuteron and the 3 Si 
scattering states. As shown by Danilov |L0[j , these components involve a pnc transition 
from a 3 Si to a 3 Pi state (the electric dipole operator does not change the spin for its 
dominant contribution) and, thus, implies aT = 0toaT=l transition. This can only 
arise from an isovector pnc force. 

The isovector pnc force is expected to be largely dominated by the pion-exchange con- 
tribution, due to the relatively low mass of the pion and, consequently, its somewhat long 



range. The expression for this force, first given in [13|, is currently written down as follows 



T r ■ QttNN h\ _ . . _ _, . 



2M N 4?rr 



(4) 



where g n NN represents the strong 7rNN coupling constant, which is taken positive here as 
usually done. Contrary to the meson-nucleon interaction, the above expression of the pnc 
potential is essentially free of any convention. The sign of V n is therefore determined by 
the sign of h\ (see Sect. 4 for more details). 

Knowing the expression of the pnc force, one can now determine the parity admixed 
components it produces in the total deuteron or 3 Si scattering states, which we write as: 

I* > d .t= Ud,t(r) + •••+< T ~^ L (* P + ^ t (r)) \S=1>, (5) 

V 2 \r p -r n \ J 

where the dots stand for the tensor component we omit for simplicity. From now on, 
particles 1 and 2 are ascribed to the proton and neutron respectively, but we nevertheless 
keep track of the isospin factor, equal to 1, to evidence the symmetry between the two 
particles. As the pnc force is a small perturbation, the expressions of the parity admixed 
components, ifjd{ r ) and ?M r )> can be obtained using the full Green's function projected 
on I = 1 angular momentum states. For the deuteron state, one gets: 

Mr) = J dr'r" G l ^(r, + W) (Mr) + ■■■), (6) 

where E = — \Ed\. A similar expression holds for the scattering state with the appropriate 
replacements. Spin and isopin dependences have been factored out so that Exp. (1) 
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(together with Exps. (2, 3, 6)) offers no ambiguity as to its sign once the sign of h\ is 
given (intrinsic phases relative to the deuteron or scattering states cancel out). 

A first but unrealistic estimate of the asymmetry, A 7 , can be obtained by using the 
following wave functions corresponding to a zero-range force together with the free Green's 
function projected on the / = 1 space: 

ltf(r) = l-^, V>?(r) = l-i $(r) = V^^, (7) 



G 



i=i 

E=0 



[r,r') = -M N — - 



(r — r 
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i=i 



(r, r') = —M 
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e~ ar (1 + ar) 



[l + ar')-e ar (1 - ar') 



2a 3 r 



/2 



r — r 




In these equations, a is related to the deuteron binding energy, a = ^Mh\E&\. The 
expression of the Green's function, only given for r' < r, allows one to precise factors. 

The calculation of integrals in Eqs. (2, 3, 6) can be performed analytically with the 
result: 



g-KNN h\M N y/2 
37r (n P - n n )(l - aa s ) 



m-n 



(niir + a) 



VVi-^ lo «( 1+ - 

a z (m 7r + a) 2cr m T 



(9) 



This expression is identical to that one obtained by Kaplan et al. when using the effective 
field theory at the leading order. In view of the extreme simplicity of the wave functions it 
corresponds to, the result can be easily improved. By using more realistic wave functions, 
one can provide insight on the corrections that should be applied to the estimate of A 7 
made by these authors. Let's notice that, due to different inputs concerning mainly the 
strong nNN coupling, the magnitude of the asymmetry given by the above equation is 
slightly larger than what they obtain. The discrepancy has not much physical relevance 
but is important in order to assess the precise validity of the approach. Thus, the value to 
be compared with more elaborate estimates presented below is — 0.187 h\ (corresponding 



to 



4tt 

23.7 fm). 



14.4, 77v = 0.7 fm" 1 , M N = 4.76 fm" 1 , fi p - fi n = 4.706, a = 0.232 fm" 



3 Size of the effect 



The size of the pnc asymmetry, Ay, has been calculated in many papers JT0|, [|, [5|, ||, 0, ||. 



They evidence discrepancies depending on whether results are corrected for the fact that 
most realistic models used in the past were not providing the right neutron-proton scat- 
tering length for the 1 Sq state. These models were assuming isospin symmetry, often 
predicting a scattering length of about — 17 fm. This represents the main source of differ- 
ence, roughly a factor 1.4 (~ )• Actually, this can be remedied by taking the value 
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of the regular transition amplitude, Ml, from the measured capture cross section, which 
offers the further advantage to also account for a 5% correction due to meson-exchange 
currents. Results also depend on whether interaction in P waves or effects due to the 
tensor force are considered. When the main ingredients mentioned above are accounted 
for ||, the uncertainty is due to the strong interaction models themselves and the precise 
way to correct for their deficiencies. Asymmetries calculated on the same footing with the 
Hamada- Johnston, Reid-soft-core and de To urreil- Sprung potential models for instance 
were found to be equal to — 0.109 h\, — 0.114 h\ and — 0.107 h\ respectively M. In all 
cases, the effect of the repulsion in the 3 Pi state (~ 25 — 30%) is partly compensated by 
an effect of about 20% due to the admixture of a 3 D\ component to the 3 Si one. Careful 
examination also indicates that part of the differences has a well determined origin: a 
slightly too large binding energy for the deuteron in the case of the Hamada- Johnston 
model and a 2% difference for the triplet scattering length for the two other models. 
The discrepancies can therefore be made smaller by using improved models of the NN 
interaction. 

The insensitivity to the strong interaction model was confirmed later on by using the 
Paris potential ||, which is a more up to date model, including an accurate treatment 
of the two pion-exchange contribution. Prior to the comparison, the result has to be 
corrected for a wrong 1 Sq scattering length however. 

The above values have a magnitude smaller than what was obtained by Kaplan et 
al., -0.187/4 with 

our inputs (see remark end of previous section). The discrepancy 
can be easily ascribed to the oversimplified wave functions, which the Kaplan et al.'s 
calculation corresponds to (Eq. 9). With respect to these ones, Hulthen wave functions 
(corresponding to a separable Yamaguchi potential) represent a significant improvement 
although they are not good ones in view of the present standard. They offer the advantage 
of providing an analytic expression for the pnc asymmetry, El, which can be used to 
estimate the effect of the finite range of the strong interaction neglected in the above 
calculation. It is given by: 



gnNN 



hlM N V2 



x 



77l,r 



+ 



r] 3 +a ijs+ijf 



(m n + a) 



+ aa t 



mi 



v a 2 (m 7r + a) 2a 3 



" ,l log(l + ^) 



m 7 



m n 



—aat 



("V + Vt) 2 2j] t {m^ + r] t f 
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4a 2 (m n + rj t ) 4a 2 (m n + i] t + a) 
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a 2 a t ( 2ml + Vt + m ^Vt _ ml ^ m^ + Vt 
rjt \ Ar)t(m w + r) t ) 2rf x 
a 2 a t f2ml - 2r] 2 + m^t - am^ + 2ar] t m ^ ~ Vt + a2 m 7T + 2r] 1 
(r) t + a) 2 \ A(rj t -a)(m n + 2r]t) 2(r] t - a) 2 m w + r] t + 

- ■ / 2ml ~ 2rj 2 + m^rjt m 2 - rj 2 ( m n + 2r], 





a 



4^(77^ + 2th) 24 ° gl + rjt > 

In the limit of a zero-range force (rjt — > 00, aat — > 1), the front factor together with the 
first line in the above equation allow one to recover Eq. (9), while the other terms vanish 
as r]^ 1 or rjl" 2 . For the value that rj is expected to take in the case of a finite range force 
(rj t = 1.386 fm -1 corresponding to a t = 5.42 fm and a = 0.232 fm -1 , and rj s = 1.203 fm -1 ), 
one gets a reduction by a factor 1.33. 
Another analytical estimate of the effect can be obtained using the Danilov's approach 
TOfl . This one, which relies on dispersion relations, only uses the knowledge of NN scat- 



tering phase shifts. It has been shown to give a good account of pnc effects in the 
neutron-proton radiative capture process, the uncertainty being of the order of 20-25% in 



absence of accidental cancellation fllB] . In this approach, the whole effect is factorized into 
two parts involving the electric dipole operator and the pnc NN force, which respectively 
are sensitive to the long and medium plus short range description of NN states. The pnc 
effect is thus determined by the strength of the pnc neutron-proton scattering amplitude, 
M N C, whose expression is given by: 

MNC = _s^KME m r-^ ( _^|)j, rW (11) 

47TV2 a t J -00 8p' b 

where 



Inserting in this equation the standard low energy effective range parametrization for 
p !_ 1 n p 

tg S t a t "'" 2 



phase shifts, r^r = — — + one gets: 



1 



F f - V) = ( for P 2 < °)> ( 13 ) 



^ aa t r t 



and from it: 



MnC ~ "lWT^r4 d9 I?TT^p (14) 

= Mn ((1 + ^x)( 2 --x + 2x 2 - 2 Ag(l + -))-- 

47rv2 aa t m n \ m w 3 x m w 

with x = 01 at r } w?r . A numerical estimate of the effect of the range is immediately obtained 
by comparing values of M^C calculated with r t = 1.765 fm and r t = 0. A reduction by a 
factor 1.6 is thus found in the case of the scattering state. This is partly reduced by the 



7 



fact that the scattering length a t is simultaneously enhanced, the real reduction factor 
being closer to 1.3. 

The apparent independence of expressions (11-14) of strong interaction models (it is 
sufficient to know the phase shifts) is partly misleading. It is likely that the off-shell effects 
that the underlying approach neglects have to be accounted for by extra contributions to 
the dispersion relation, Eq. (11). 

Estimates due to a better description of the S states given above do not take into 
account the repulsion at short distances. From the comparison with results from realistic 
NN interaction models, an extra reduction by a factor ~ 1.15 — 1.20 is found for this part. 
Some reduction by a factor ~ 1.25 — 1.30 is also due to the repulsion in the 3 Pi-state. 
Thus the total reduction of Ay with respect to the value — 0.187 h\ amounts to a factor 
~ 1.9 — 2.0. The size of the effect is larger than the 30% referred to in ref. [[H]]- With this 
respect, we would like to mention that this uncertainty is almost reached by the effect of 
neglecting the finite range of the interaction. The value of the product aa t , which is set 
to 1 in the zero-range limit, is actually larger by a factor 1.26. The deuteron asymptotic 
normalization, \^2a, is found smaller than the experimental one, 0.884 fm~a by a factor 
1.30. As these discrepancies concern long range properties, one should not be surprised 
by the significantly larger discrepancy we found for the pnc effect. This one depends on 
the short range properties that are much less well described by the present use of effective 
field-theory methods. 

The agreement obtained for the total capture cross section was considered by Kaplan 



et al. as a support for their approach [ll|. In view of the previous comments, it looks 
somewhat accidental. It results from a cancellation of two effects: a long range contri- 
bution that is too small and a shorter range contribution (further enhanced for the pnc 
case) that is too large. 



4 Sign of the effect 

Signs in the field of parity violation are an important information but have not always 
received the attention they deserve. The size of the measured pnc asymmetry in pp 
scattering at low energy has been reproduced in many papers for instance. A careful 



examination nevertheless showed that ten of them were actually missing the sign [17] . In 
view of this, we here precise ingredients that enter into the determination of the sign of 
the pnc asymmetry, A 7 , together with its expression. They concern the pnc potential, V^, 
the electromagnetic transition operators, and the calculation of the asymmetry itself. 
The determination of the sign of the pion-exchange pnc force, Eq. (4), has been made 



for the first time in ref. JT2|. It makes sense as far as the expression of Pauli matrices 



and the relation p = —iS7 are universally accepted. The sign is then related to the sign 
of the product of the strong and weak 7rNN couplings. Individually, the sign of each of 
them depends on many conventions: intrinsic phase of the meson field, definition from an 
Hamiltonian or from a Lagrangian, expression of the 7 matrices, especially the 75 matrix 
which makes the difference between a pc and a pnc interaction, spin-| fields referring to 
particles or anti-particles, metric, ... While it is possible to get rid of all these conventions 
with some care, this is not sufficient to determine the sign of the product of the couplings. 
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The last step is achieved by noticing that the calculation of the pnc coupling, h\ y using 
PCAC and current algebra, is involving a factor proportional to the strong coupling, hence 
a dependence of (j-knn h\ on the quantity gl NN , which is always positive [Q. The sign 
of the three types of contributions considered in DDH (sum rule, "parity admixture" into 
the wave function and factorization; see Fig. lc, lb, la in the Appendix) could thus be 
determined. Taking as positive the strong coupling appearing in the pnc potential, Eq. 
(4), as done in DDH and subsequent papers 0, ||, the estimated weak coupling h\ and 
its range given in these works correspond to a positive value. 

The electric El and magnetic Ml operators are a possible source of sign mistake. They 
can be derived from the electromagnetic Hamiltonian, which may be written as: 



{Pi+Pf) 



1 + T Z 



2M 



A' 



■2M 



+ 



T 



N 



2M N 



)j Oj . q x e + 



(15) 

where the dots account for meson-exchange currents required to ensure gauge invariance. 
The photon polarization is represented by the vector e and its momentum (energy) by 
q{q°). An alternative expression, which is of more direct usefulness for the present work, 
relies on the Siegert theorem to account for the dominant contributions due to exchange 
currents (represented by dots below), 



p . e 



.1 +T Z 



2M 



N 



), + .... 
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.1 + T Z 



It is given by: 



H, 



).e±l~ rrr (CT„ 



2 2M 



N 



o n ) . q x e 



+ .... 



(16) 



where the dots, here and below, represent contributions that are irrelevant for our purpose. 
Different signs appear in Eq. (16) as well as Eq. (15), depending on whether an absorption 
or an emission of the photon is considered (or on whether one calculates the matrix element 
< i\H\f > or < f\H\i >). The present ones are given for an emission (radiative capture) 
and respectively correspond to < i\H\f > (< np\H [d'y >) and < f\H\i > (< dj\H\np >) 
for the upper and lower sign. Either way leads to the same result, but we nevertheless 
emphasize the point as the electric dipole operator, r, in Eq. (16) has not always been 
used with the appropriate sign for the energy factor q°, giving rise to different signs for 
the parity violating asymmetries || . 

Using the expression of the electromagnetic interaction, Eq. (16), the expression of the 
total wave function, Eq. (5), together with Eq. (6) for the parity admixed component, 
one obtains the transition matrix operator in terms of the matrix elements, Ml and El, 
previously defined, Eqs. (2, 3): 



T(np -> di) oc (±i Ml (o p - a n ) .qxe + El (a p + a n ) . (q° e- qe°) + ... 



(17) 



The last term has been completed with a contribution from the charge density so that to 
emphasize its gauge invariance. Using the relation, i (dp — dn) = (dp x d* n ) P a , one obtains 
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an alternative expression. Less dependent on conventions as far as the relative sign of Ml 
and El contributions is concerned, it is given by: 

T(np -f dj) oc (Ml (dp x a n ).qx e + El (dp + a n ) . (q° e - qe°) + ....) . (18) 

The angular distribution of the emitted photons with respect to the neutron polarization 
n is obtained by taking the matrix element of the transition operator, T, corresponding to 
a given neutron polarization along this direction, multiplying it by its complex conjugate 
and summing over all polarizations in the final states. This amounts to calculate: 

a(6) oc ETr. ( (1±^) T (5 + ^) T+ ) 

oc (Ml 2 -2 Ml El n.q + ...), (19) 
from which Eq. (1) immediately follows. 



Present estimates and constraints for the pnc cou- 
pling, hi 



There are in the literature numerous papers dealing with the estimate of the pnc coupling, 



h\, in the standard electro- weak interaction model 



12, 19, 14, 20, 21, E 



26 1 . The most comprehensive work is still probably the DDH one \ T% . Contributions 



considered since then most often correspond to one of those estimated in this earlier work. 
While the improvement is not always obvious, at least, they provide a better estimate of 
the uncertainty that should be ascribed to each of them. Due to expected cancellations 
between different contributions, picking up one of them is however misleading. Some 
details and a discussion are given in ref. |§ as well as in the Appendix. One of the 
features evidenced by the results, summarized in DDH, was their strong dependence on 
the factor, K, that accounts for strong interaction effects. Since then, due to improvements 
in determining the strong QCD coupling, the value of K has decreased from K = 6 to 
K = 3 for an energy scale of /x = 1 GeV. It was also realized that there was some 
cancellation between the effect due to a dependence of the K factor on this (arbitrary) 
energy scale and that one for the quark masses p7| |, with the result of making estimates 
of pnc meson-nucleon couplings less dependent on it, as it should be. The introduction 
of a contribution implying strange quarks in the nucleon could have the same effect for 
another contribution involving strange quarks in the hadronic weak interaction. Thus, an 
up-dated range for the coupling would now be: < h\ < 2.5 1CT 7 (see [^] for an analysis 
of some uncertainties). Due to large cancellations in the total contribution involving u 
and d quarks (see Table 4 of ref. || and Appendix), most of the contribution arises from 
the part of the weak interaction involving strange quarks. It is interesting to notice that a 
recent estimate of the same contribution |26| falls in the above range. Other estimates fall 
outside. While that one for the strange quark in ref. |25[] is not a serious source of concern 
in view of its very rough character, those based on QCD sum rules pp| , |23| , |2~4] , which in 
any case ignore strange quarks, are. It is not clear which contribution they correspond 
to in the DDH picture. Let's finally mention that all studies give positive values for h\. 



10 



From the examination of the papers however, we did not get the certitude that the sign 
was carefully determined (assuming it was looked at). 

Phenomenological constraints on the coupling, h\, have been recently discussed in the 
literature [^, 59]. To explain the recent measurement of the anapole moment in 133 Cs 
30fl , it has been suggested that a large value of the coupling would be required |31j, |32|] , 



exceeding the upper limit determined from the experiment performed at ILL in the past 
|| ( h\ = (1.4 ± 4.4) 10~ 7 , in absence of other contribution). An isoscalar pnc force 
enhanced in the nuclear medium could do as well with a coupling that would be smaller 
and essentially in agreement with what the absence of observed effect in 18 F indicates, 
h\ < 1.3 1CT 7 [Q. It is remembered that the calculation of the nuclear pnc matrix element 
in this case can be checked by using information from first forbidden (] decay from 18 Ne. 
Although it has not the same probing strength as in 18 F, the absence of effect in 21 Ne 
and in 93 Tc also points to a small coupling if one discards it results from accidental 
cancellations [[3], |33| . Let's remind that the present theoretical contribution of the pion- 
exchange to the circular polarization of photons emitted in the transition, ^~ — > ^- in 
21 Ne is expected to be about \P*\ = 10 5 \ h\\, whereas the measurement puts the limit 
P 7 = (0.8 ± 1.4) x 1(T 3 . 



6 Conclusion 

We have shown that, far from representing a progress, the recent calculation of the pnc 
asymmetry in n+p —>■ d+7, based on using effective field-theory methods at leading order, 
is a step back to pre-historic times of the field. While this approach looks appealing at 
first sight, it evidences large discrepancies with previous calculations. Realizing that these 
latter rely on much more elaborate descriptions of the deuteron and scattering states, they 
can be used to assess the validity of the new approach as used until now. 

From the comparison we made, the uncertainty of 30% ascribed to it by the authors is 
reasonable for the calculation of observables that have a long range, such as magnetic or 
electric transitions. In the present case however, part of the effect involves a weak pnc 
interaction which has a shorter range. Examination of their expression for the asymmetry 
Ay, which corresponds to an attractive zero-range force, shows that it behaves like — , 
while for reasonable potentials one rather expects -\ (all other quantities being supposed 
to be small). This part of the calculation is overestimated by a factor of about 2. The 
discrepancy has its origin in the treatment of nucleon-nucleon correlations and especially 
their well known short range repulsive part, which prevents two nucleons to have an 
infinite probability to be close to each other (phase space factor put apart) as the present 
calculation by Kaplan et al. supposes. It decomposes into a factor roughly given by 
(1 + c - jp* ) ~ 1.5 — 1.6, representing the combined effect of a finite range attraction 
together with a short range repulsion for the force in the 3 Si state (r t represents the 
effective range and c ~ 1.3) and a factor ~ 1.25 — 1.30 representing the effect of long and 
short range repulsion in the 3 Pi states. The overestimate is actually smaller by a factor 
~ 1.15 due to a contribution from the tensor force that is constructive in the present case, 
itself reduced by the correction of mesonic exchange currents to the regular transition, 
Ml. While the new result by Kaplan et al. lets hope that the pnc asymmetry, Ay, could 
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be quite large in view of the error ascribed to the theoretical method, it thus appears that: 
1) the error with respect to the actual value is underestimated, 2) this value is slightly 
outside the range they provided, in the lower part. 

The underestimate of the error indicates that the parameter which, in the approach, 
governs the expansion at different orders does not provide a good magnitude of their 
relative contributions. This conclusion is supported by the omission of the dominant 
contribution to the deuteron anapole moment in |34| (corrected in |35| ) which suggests 
that the determination of what is relevant at some order is not so straightforward. This 
does not preclude the use of effective approaches for describing pnc effects however. Such 



approaches were initiated in the field by Danilov with a different foundation |T(J and 
extented to nuclei in || . They can give a good account of pnc effects in a large variety of 
low energy processes. 

The determination of the sign of the asymmetry requires some care, especially to get rid 
of numerous conventions that appear at different steps of the calculation. Assuming that 



the authors of ref. JTT[] have relied on commonly used definitions for writting down the 
expression of the pion-nucleon coupling, we found that the pion-exchange potential which 
it gives rise to has a sign opposite to the current one, thus producing an asymmetry with 
the sign different from that one generally referred to. After correcting their result for 
this different choice, the determination of the sign by these authors confirms the earlier 
determination, thus providing the independent check that they were calling for. 

Concluding, there is a long way ahead before the effective field-theory approach to the 
calculation of the asymmetry A 1 be able to compete with more elaborate calculations. For 
the pion-exchange contribution, these ones converge to the value = —(0.11 ± 0.01) h\. 
The sign refers to the current expression of the 7r-exchange pnc potential and the error, 
estimated largely, results from the comparison of the most realistic descriptions of the 
process. 



Appendix: More about h\ 



Details on the theoretical estimates of the pnc coupling constant, h\, are given here. As 
well known, this coupling depends on the isovector part of the hadronic weak interactions, 
which reads in the "free" case: 

H pnc (AT = l) = -j= (^sin 2 6 c [rfuvrfffcs + s^ l5 uu llx s - u -> d] 

-^sin 2 6 w sy^s (uj^u - dy^d) 



1 - 2sin 2 9 



-^sin 2 6 w (uYu + d^d) (uj^u - dj^ 5 d) \ . (20) 
In practice, it can be calculated from an effective interaction which incorporates effects 
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Figure 1: Quark processes contributing to the pnc meson- nucleon couplings. Diagrams 
labelled a, b, c, d, arise respectively from the factorization approximation, parity admix- 
ture into the nucleon, first order and part of the second order sea quarks. Diagram c) 
has been corrected with respect to the one presented in Ref. [[|. In all diagrams, the 
box represents the pnc quark interaction. Its effective character can be seen from the 
examination of the last two diagrams. By extending the box, it is possible to have it 
to incorporate that part involving gluon exchange so that to get a diagram a) type. On 
the other hand, the qq pair, which the gluon couples to, contains in particular a ss pair, 
colored in c), partly uncolored in d). This last diagram involves a sub-diagram similar to 
that one responsible for reducing the contribution of quarks to the nucleon spin. 
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K = l + ^L^H^f)- (21) 



due to strong interactions (QCD). In the four-flavor case, this interaction involves the 
factor, 

2 5g 2 (/i) M 2 ,, 
48vr 2 1 /i 2 

The quantity, </ 2 (/i), is the QCD running coupling constant at the energy scale /i. This 
energy scale is in principle arbitrary, but an optimal value often quoted in the literature 
is fi = 1 GeV. The determination of this coupling has varied (and improved) with time. 
It has dropped from the value 1.0 used in DDH, corresponding to K = 6, to a value 0.4, 
corresponding to K = 3. For this value, the detail of the various contributions to the 
coupling constant h\ estimated in DDH reads: 



= 1 + 3.8 - 3.3 + 0.5 + 3.3 ) = 2.0 10~ 7 . 

ss + ss + {uu + dd) + par. admixt. + fact. (22) 

The sign in this expression refers to a positive sign of the strong coupling, g^NN, ap- 
pearing in Eq. 4. Otherwise, all the contributions are normalized to the contribution 
of the charged current, /° = 0.38 10~ 7 , which was calculated first in [H| and involves 
strange quarks (this is remembered below the corresponding contribution). The other 
contributions in Eq. (p2[) arise from the neutral current. The second and third contribu- 
tions, like the first one involves sea quarks, respectively strange and non-strange. They 
can all be traced back to the diagram c of the figure. The fourth contribution involves 
"parity- non-conservation" in the nucleon wave function (see diagrams b of the figure). 
The last contribution corresponds to the factorization approximation (diagram a of the 
figure). 

Concerning the order of magnitude, the first term in Eq. Q2"2] ) contains a factor sin 2 9 C ~ 
0.05, explaining for some part its small contribution. It arises from the first line in Eq. 
(20). The second one essentially involves a factor 1 — 2sin 2 6'yi/ ~ 0.5, which makes its 
contribution larger than the first one. It arises from the second and third lines in Eq. (20). 
Being of the second order in gluon exchange, instead of the first order, the enhancement 
is not as large as one could naively expect from the comparison of the factors mentioned 
above. The three other terms in Eq. ( 22|) contain a factor sm 3 w ~ 0.08. They only 
involve u and d quarks and arise from the fourth line in Eq. (20). While examination of 
factors give some insight on the magnitude of each term, making a precise prediction is 
difficult due to the presence of a destructive contribution. Depending on the prejudice on 
the size of the various contributions, one can reasonably expect the coupling constant, h\ 
to be in the following range, (—1, 3) x 10~ 7 . 

It is instructive to make a detailed comparison of the above results with other works. 
Prior to DDH, Weinberg considered the contribution from neutral currents involving 



strange quarks fll9| . The value he got was much larger than obtained above but was 
calculated in the limit K — > oo. The finite value of K reduces his estimate from about 
12 to 3.8 in Eq. (p2[). A contribution involving sea- but non-strange quarks was also 



considered by Gari and Reid [18]. After correction for a factor -4, their contribution is 



decreased to provide a value that is accounted for by the negative contribution -3.3 in 



Eq. (P2j). Dubovik and Zenkin [21] considered contributions that can be compared (and 
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do compare) to the fourth and fifth numbers in this equation. Using QCD sum rules, 



Khatsimovskii |20| made an estimate which turns out to be larger than in Eq. (|22|) and 
close to the "best guess" DDH value. In principle, it should involve some of the three 
last contributions in this equation (see comment below). Kaiser and Meissner |[22| , using 
a quite different approach based on a chiral Lagrangian, calculated a contribution that 
involves non-strange quarks and got a small number for the coupling h\ . This one com- 
pares reasonably well with the sum of the three last contributions in Eq. (|22|) . With this 
respect the possibility that one of the contributions has a destructive character is impor- 
tant and should deserve some attention. The role of strange quarks, neglected in the last 
four works just mentioned above but taken into account in earlier works was reminded 



by Kaplan and Savage ||25|| . Their estimate agrees with that one made by Weinberg and, 
in view of its very rough character, can be considered as consistent with the contribution 
proportional to 3.8 in Eq. (^). Another QCD sum rule estimate by Henley, Hwang and 
Kisslinger |23) , always neglecting strange quarks, first turned out to be in agreement with 
what was expected from the three last numbers in Eq. (E3). After correction for the 



relative sign of two contributions (see errata in |f23] and p4j1), this estimate came closer 



to that one made previously by Khatsimovskii. More recently, Meissner and Weigel [BHj 
considered the contribution of strange quarks in a chiral model. Their estimate is close 
to that one given by the number 3.8 in Eq. (|22l). 

While most contributions considered in the literature have a counterpart in the DDH 
framework, there is a serious problem with QCD sum rule based calculations. One can 
imagine that they take into account a new contribution, but one would like to know 
which one in terms of quarks, knowing that DDH accounts for the simplest ones depicted 
in Fig. 1 (a,b,c). The other possibility is that something is missed by these approaches. 
With this respect, we would like to mention that a recent examination of the contribution 
involving sea quarks (Fig. la) led us to a contribution with a sign opposite to what is 
suggested by the analysis of hyperon decay amplitudes, used in DDH. In such a case, the 
two contributions due to non-strange quarks involving the qq vacuum expectation value, 
< g.s.\ qq \g.s. >, (3rd and 5th terms in Eq. (|22"D), would add together. However, it seems 
that such an estimate misses a contribution that is required to ensure the conservation of 
an axial current coupled to the gluon field (there is some similarity with the calculation 
of an anapole moment). Enforcing this conservation might lead to a contribution with an 
opposite sign, in agreement with Eq. (p^). 

As mentioned above, the energy scale, which the K factor depends on, is arbitrary. 
Two contributions in Eq. fl2"2| ) evidence a strong dependence on this factor K and, through 
it, on the energy scale (see U). This one, which determines the frontier between strong 
interaction effects incorporated separately in the effective weak interaction and in the 
description of hadrons, should not intervene in the determination of physical quantities. 
In the case of the fifth contribution to Eq. (^) (factorization), this result is qualitatively 
achieved by cutting off the high momentum part of the effective interaction between 
quarks (as is expected), which, most probably, is equivalent to the effect of introducing 
a dependence of quark masses on the energy scale [071] . In the case of the strange quark 
contribution (second term in Eq. (f22|)), which vanishes in absence of strong interaction 
(K = 1), the energy scale independence requires to introduce a strange quark-antiquark 
pair, ss, in the description of hadrons. A contribution that may play some role in this 
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respect is represented in Fig. Id. It has some relationship to the contribution that tends 
to reduce that part of the nucleon spin carried by quarks. 

Concerning notations, let's mention that the coupling h\ used here is denoted h\ NN) 
or f n in other works. The notation /„. is confusing as it represents another quantity 
frequently used in the field of strong interactions. The coupling, h^, is also related to the 
quantity introduced in pj, Hi = = 2.38 /£. 
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